On the spectral gap in Andreev graphs 

Holger Flechsig and Sven Gnutzmann 

Abstract. We introduce Andreev scattering (electron- hole con- 
version at an interface of a normal conductor to a superconductor) 
at the outer vertices of a quantum star graph and examine its ef- 
fect on the spectrum. More specifically we show that the density 
of states in Andreev graphs is suppressed near the Fermi energy 
where a spectral gap may occur. The size and existence of such 
a gap depends on one side on the Andreev scattering amplitudes 
and, on the other side, on the properties of the electron-electron 
scattering matrix at the central vertex. We also show that the 
bond length fluctuations have a minor effect on the spectrum near 
the Fermi energy. 



1. Introduction and physical background 

The spectra of superconducting-normalconducting (SN) hybrid de- 
vices such as Andreev billiards has been a topic in physics for some 
time (see the recent review pLj for an overview and further references). 
An Andreev billiard consists of a finite normalconducting region with 
a boundary that is partly an interface to a superconductor. An elec- 
tron with energy near the Fermi energy cannot propagate into the 
superconductor and it is trapped in a finite region similarly to a usual 
quantum billiard. However, the interface to the superconductor leads 
to an effective coupling of electron and hole dynamics through Andreev 
scattering. If an electron hits the interface to the superconductor it is 
reflected back into the billiard as a hole in the opposite direction of the 
incoming electron (this has been called "retroflection" to distinguish it 
from the specular reflection at an interface to an insulator) . Physically 
the electron enters the superconductor where it excites an electron-hole 
pair. It then builds a Cooper pair with the electron while the hole is 
ejected back into the billiard. Quasiparticle excitations in such a device 
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are described by the Bogohubov-de Gennes equation 

/-V2 + V{x) - /i A(x) \ 

(1-1) , w2 T.. ^ U(x)=i?V(x). 

\^ A(x) V^-V{x.)+fiJ 

Here, V{x.) is a potential, fi = kp is the Fermi energy [kp is the Fermi 
wave number), and E the energy measured from the Fermi leveL The 

wave function i/'(x) = ( ^^^j 1 has two components: u(x) is the elec- 
tron component and f (x) the hole component. The electron and hole 
components are coupled by the pair potential A(x). The pair po- 
tential vanishes in a normalconducting region. In superconducting 
regions A(x) 7^ and waves can only propagate for E > |A(x)|. 
In (11. ip we have assumed time-reversal invariance which means that 
there is no magnetic field and the pair potential is a real function. 
/m(x)\ 

If i/'lx) = ( ^^x) ] ^^^^^^ Bogoliubov-de Gennes equation at en- 

ergy E then Ci/'(x) = ( 1 is a solution to energy —E. This 

shows that the spectrum is symmetric around E = (that is around 
the Fermi energy) and reveals that these systems belong to one of 
the seven non-standard symmetry classes [2] (the three standard or 
Wigner-Dyson symmetry classes classify quantum systems according to 
their behaviour under time- reversal and spin rotation). The negative 
spectrum describes hole-like quasiparticle excitations (while the spec- 
trum is negative, the physical energy of such an excitation is positive, 
of course) and the positive spectrum describes electron-like quasiparti- 
cle excitations. 

It has been known for some time that Andreev reflections reduce the 
density of states near the Fermi energy and various universality classes 
have been defined [l[3l[4l[5l[6l[7lll[9l[lQ]and related to quantum 
chaos. Spectral gaps have been found in irregularly shaped Andreev 
billiards for which the classical dynamics of the normal billiard (with 
specular reflection at the superconducting interface) is chaotic. Though 
Andreev scattering destroys chaos the Lyapunov exponent can be felt if 
the mean time between two subsequent Andreev reflections is large 
- the gap that was found is of the order h/tA (in general much larger 
than the mean level spacing) [31, |4l, [5] . In contrast an integrable nor- 
mal billiard only leads to a suppressed density of states near the Fermi 
energy [T], [8] . Other universality classes such as magnetic billiards may 
have a chaotic dynamics even in the presence of Andreev reflections - 
in that case the effect of Andreev scattering is universally described by 
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appropriate Gaussian random- matrix ensembles [2], 110] . However the 
effect of Andreev scattering on the spectrum is less strong and decays 
after a few mean level spacings away from Fermi energy (which can 
only be seen after sufficiently averaging over system parameters). 
Andreev scattering has been introduced to quantum maps [TTj and 
quantum graphs |12] . The latter work investigated a class of Andreev 
graphs that corresponds to chaotic combined electron-hole dynamics 
with universal spectral statistics for non-standard symmetry classes. 
Chaos is however suppressed in a non-magnetic Andreev billiard. Here 
we will discuss the effect of Andreev reflections in a different class of 
graphs that is more closely related to non-magnetic Andreev billiards. 

In Section [2] we will start our discussion with a simple model: a 
single lead with one interface to a superconductor at one end. This will 
serve as a motivation for our construction of Andreev graphs in Section 
13. 1[ In the rest of Section [3] we will discuss the spectra of various 
ensembles of Andreev graphs - partly analytically, and numerically 
where analytical results could not be obtained straight forwardly. 



2. A hybrid SN lead 

Most of the concepts we will use later are most easily introduced 
for a one dimensional electronic hybrid SN lead. 



2.1. Scattering from a SN interface. Let us start with dis- 
cussing scattering from a SN interface at x = 0. Let us assume that 
the pair potential is constant on the superconducting half line x > 
where A(x) = Aq while A(x) =0 in the normal conducting half line 
X < 0. At the interface we include a potential barrier V{x) = VqS^x). 
The scattering problem can be solved exactly for any set of parameters. 
We are interested in the case < Aq in which there is no propagation 
in the superconducting half line. On the normalconducting half line 
a scattering solution of the Bogoliubov-de Gennes equation takes the 
form 



(2.1) 



NC'. 



—ikf^x 



'-''e,out n — 



+ 



V '"e 



where = \/fx + E and = yJi — E are the wave numbers of elec- 
trons and holes. Note that a hole propagating to the right (incoming 
wave) is described by e"*'^''^. In (12. ip the outgoing coefficients ae//i,out 
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are related to the incoming coefficients by the scattering matrix S{E) 

«e,out \ / See{E) Seh{E)\ f a^^i^ 

ah ,out / \SheiE) ShhiE) I \a/jii- 



(2.2) 



The scattering matrix can be obtained by matching the solution in the 
normalconducting part with a decaying solution in the superconducting 
part at X = 0. The regime E Aq <ti ^ = kp will be relevant to us and 
simplifies the expressions considerably. Let us formalise this regime by 
using the observation that the difference of the electron wave number 
and the Fermi wave number is proportional to the energy 
■\/kp + E—kp = ^+0{E^/kp). Keeping this difference fixed as fci? 
oo defines an interesting asymptotics. We thus introduce k = ^ as the 
new (rescaled) energy variable. All other system parameters have to be 
rescaled appropriately. For the pair potential we write Aq = k]^'^5^^°' 
where e < a < 1 — e which ensures E <^ Aq fi. For the potential 
barrier Vq = kpVo is the only way to scale the strength such that the 
resulting limiting scattering problem still contains a finite barrier. The 
leading order 5*^ of the scattering matrix as kp ^ oo (keeping n, 6 and 
Vq fixed) only depends on the potential barrier strength vq, explicitly 



(2.3) Sa 
where 



(2.4) t = 2 e 



it -Vl^e-'f^ 

Vq + 2i 



ifS 



2 + 



If there is no potential barrier vq = the scattering matrix Sa_ de- 
scribes a pure Andreev refiection where an incoming electron (hole) is 
refiected as a hole (electron) and a phase factor She —i is acquired. 
In the opposite limit of a large potential barrier f o ^ oo one has pure 
electron-electron (hole-hole) scattering with See, Shh ^ — 1 which is 
equivalent to Dirichlet boundary conditions. 

The leading asymptotics for the wave function (12.11) is obtained by re- 



placing replacing S ^ Sa and Note that this 



contains fast oscillations oc e^^''^. In the following we will always as- 
sume that kp is finite but large enough that it is justified to consider 
the leading asymptotic term only. 

2.2. Excitation spectrum of a finite lead with one interface 
to a superconductor. We can now discuss the discrete excitation 
spectrum {k„ = ^} of a finite normalconducting lead of length L 
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(— L < X < 0) which is attached to a superconductor at x = and to 
an insulator at the other end x = —L. We assume Dirichlet boundary 
conditions i/?nc(^ = ~L) = at the interface to the insulator. Let us 
define the quantum map 

(2.5) U{k) = SaT{k)SdT{k) 

where 
(2.6) 

= -ij ^"""^ ^^'^^ = e-^C^^--^)^ 

Here T{k) contains the phase factors acquired by the electron and 
hole components when propagated along the lead from one side to the 
other. The quantum map contains four steps: i. propagation from 
the superconducting interface to the insulating interface with T{k), 
a. scattering from the insulator with So, Hi- propagation back to the 
superconducting interface, and iv. scattering from the superconducting 
interface with Sa- If the quantum map U{k,) has an eigenvalue unity 
then K belongs to the excitation spectrum. The quantum map for 
the lead is periodic U{k + = U{k) - the spectrum has the same 
periodicity and it suffices to know the first two positive eigenvalues ki 
and K2. Explicitly, the quantum map is given by U{k,) = c^^^'^Uq with 



(^•'^) ^0 — \ ^i2kFL ^1 _ f2^-i2kpL+/3) 

and has the eigenvalues 

(2.8) Ml = e*(2LK-0) ^j^j ^ ^i{2L.+<p) 

where 

(2.9) = arccos (cos{2LkF + (3)Vl-t^^ G [0, tt] . 

The first two positive eigenvalues in the spectrum of the lead are ki 

2L '^2 - 



^ and K2 = ^ We arrive at the corresponding density of states 



n=— oo ^ ' ^ ^ 

We can now consider the mean density of states {d{K,))kp averaged over 
different values of the Fermi wavelength kp where (0)^^ = limx^oo Iko~^^ © ^^f- 
The result will only depend on the barrier strength vq or, equivalently, 
on the probability pa = to be Andreev scattered at the supercon- 
ducting interface. From (12.91) one observes that is always confined to 
the interval arcsin t < cj) < -k — arcsin t and every value in that interval 
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Figure 1. Averaged spectral counting function for a 
single lead (see text) for some values of the probabil- 
ity of Andreev scattering p^. The same plots also apply 
to homogeneous Andreev star graphs with B bonds and 
negligible bond length fluctuations (see Section [3.2.11) . 



is assumed for some kp. As a consequence the mean density of states 
contains gaps of size 2Ag where 



(2.11) 



a.t K = nn {n = 0, ±1, ±2, . . . ). The average over the Fermi wavelength 
kp can be performed straight forwardly by replacing kpL + (3^9 and 
averaging over one period of 9 



(2.12) 



{d{K)), 




sin(2LK) 



for < K < ™^ 

arcsin t ^ ^ ^ tt— arcsin t 



2L 



2L 



for * < < ^ 
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Figure [H shows the corresponding averaged spectral counting function 
(2.13) 

for < K < S£|hU 
iarccos { ^^7^) for ™^ <k< 




r TT— arcsin t ^ ^- ^ tt 



for various values of the probability pa = t for Andreev scattering 



3. The spectral gap in Andreev graphs 

In this chapter we will generalise the discussion of the lead with 
one superconducting/insulating interface at the ends to Andreev star 
graphs. A star graph Qb consists of one central vertex and B outer (or 
peripheral) vertices which are connected to the centre by B bonds (or 
edges). We assign a length Li, {b = 1, . . . , B) and a position variable 
Xb E [0, Lh] to each bond. The star graph can then be quantised (and 
become a quantum star graph) by defining a Schrodinger operator on 
the graph (see |13L I14j and references therein) . An Andreev star graph 
can similarly be constructed by defining a Bogoliubov-de Gennes op- 
erator on the graph. 

Alternatively quantum graphs can be defined using the scattering for- 
malism of Kottos and Smilansky [15] . Their approach has been used 
to construct Andreev graphs [12] (and other classes of quantum graphs 
that contain particle- hole symmetries). We will follow the latter ap- 
proach and assume that the peripheral vertices of the star graphs have 
an interface to a superconductor such that Andreev reflections can oc- 
cur with a certain probability. The rest of the graph (the bonds and 
the central vertex) are normalconducting such that electron and hole 
components are not coupled. The construction takes care that the 

symmetry i/'(x) = ^ Ci/?(x) = ^ ^^x)* ^ ^^^^"^^^ transforms 

positive energy states to negative energy states is also obeyed by the 
Andreev graph. Indeed the lead that we have discussed in the previ- 
ous section is already an example of such an Andreev star graph with 
-8 = 1. We will be interested what happens to the gap for some 
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ensembles of Andreev star graphs and how the properties of the den- 
sity of states depend on the properties of the central vertex and the 
Andreev scattering amplitudes at superconducting vertices. 

3.1. Construction of Andreev star graphs. Andreev star graphs 
can be quantised using a quantum map IA{k) in a similar way as we 
have treated the lead in the previous section. For a star graph with B 
bonds the quantum map is a 2B x 2B matrix of the form 

(3.1) U{k) = SaT{k)ScT{k) 

where 



(3.2) ^('^)-\^ e-*^^^-'^)^ 

C =diag(Li, . . .,Lb) 

is a diagonal matrix which describes the propagation of an electron 
(first B entries) or hole (last B entries) from one end of each bond to 
the other {C is the diagonal B x B matrix which contains the lengths 
of the bonds), Sa describes the scattering at the peripheral vertices 
and Sc the scattering from the central vertex. The (partial) Andreev 
scattering at the peripheral vertices has four blocks 



(3.3) S, 



-r —it 
-it —r 



which contain the two diagonal matrices 

t =diag(ti, ...,tB) 

(3-4) / / 

r =diag(-^l-t?,...,^l-t|) . 

The real parameter > is the amplitude for Andreev scattering at 
the 6-th vertex. We will not allow negative values for the Andreev 
amplitude her^. Comparing this to the scattering matrix at a super- 
conducting interface in the previous section one may observe that we 
omitted the phase factor e^^^ in the electron-electron/hole- hole ampli- 
tude. These turn out to be irrelevant once we average the resulting 
spectra over the Fermi wave length (in the same way as they have been 
irrelevant for the lead - here we omit them from the start). 
The bonds and the central vertex are assumed to be normalconduct- 
ing. The electron-hole symmetry of the Bogoliubov-de Gennes equation 



^Andreev graphs with mixed signs for the Andreev scattering amphtude belong 
to a different universaHty class |12| - in an Andreev billiard this would correspond 
to having two interfaces to different superconductors whose pair potentials have a 
relative phase e^^ = — 1. No strong spectral gaps can be seen in such a regime. 
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and time-reversal invariance then restrict the scattering matrix of the 
central vertex to the form |12] 

a 
a* 



(3.5) 5, 



c 



where a is a symmetric unitary B x B matrix. 

We will write Qb (W('^)) for an Andreev star graph with B bonds and 
quantum map 1J{k). The spectrum of the Andreev star graph consists 
of the discrete (positive and negative) values n = Kn for which the 
quantum map has an eigenvalue unity. Equivalently, the spectrum is 
given by the zeros of the characteristic equation 

(3.6) spec [GBiUiK))] = {kER: ^g^(ui.))i^^, 1) = O} 
where 

(3.7) ^efl{w(«))('«, A) = det (1 - A 

The spectrum (13.61) is symmetric around k = - if /€„ is in the spectrum 
so is — This can be seen from the electron-hole symmetry 

(3.8) u{Kr = (^\ ~0 

It will be convenient to define the density of states as 

(3.9) diK) = i^5(«:-«:„) 

n 

which differs from the conventional form by the factor 1/B. Accord- 
ingly we define the spectral counting function 

(3.10) N{k)= I d{K')dK' = ^^Q{k- Kn) 

(where we have assumed that k = is not in the spectrum and 0(x) 
is the Heaviside step function). The Kottos-Smilansky trace formula 
[15] expresses the spectral counting function (and hence the density of 
states) in terms of the characteristic equation 

N{k) =Arweyl(/€) + A^osc(«;) 

r ^ 2tr(/:) 

(3.11) N^ey\{K) =-g^'^ 

A^osc(ft:) = ^ lim Imln^g^(2^(«))(/t, A) . 

The contribution Nosc{i^) is an oscillating function of the spectral pa- 
rameter K while iVweyi('^) describes the mean increase of the number of 
states over a large interval of the (rescaled) energy k. The oscillating 



10 



HOLGER FLECHSIG AND SVEN GNUTZMANN 



part of the counting function does not necessarily vanish if averaged 
over some system parameter at fixed k. 

3.2. Andreev star graphs with negUgible bond length fluc- 
tuations. Let us introduce the bond length fluctuations by writing 

(3.12) Lb = L + 6 Lb where ^^6^ = 

b 

where L is the mean bond length. If the fluctuations 5Lf, are very 
small, and for small enough k one has e'^^'''^ ~ e^^'^. Indeed we will 
only be interested in the spectrum on a scale \k,\ < C-^ (such that 
^WeyK/^) ^ C and C is of order unity) and assuming SL{,j^ <^ 1/C 
we will set e^^^^'^ i— > 1. As a consequence U{k + f ) = U{k) is periodic 
(and hence is the spectrum). Also, the spectrum and eigenfunctions of 
the Andrccv graph can be constructed straight forwardly from the 2B 
eigcnvahics and eigenvectors of the quantum map W(/t) at k = 0. 
The small fluctuations in the bond lengths can however not be neglected 
in the phase factors containing the Fermi wavenumber kp - indeed we 
want to average spectral functions over arbitrary large values of the 
Fermi wavenumber. A further significant simplification is obtained by 
assuming that all bond length are rationally independent. By setting 
^ikpLi, gi^j j^g^y then replace an average over the Fermi wavelength 
by an average over a S-torus 

Writing 6 — diag(^i, . . . , 9b) we arrive at the following class of ensem- 
bles of quantum maps for an Andreev star graph with B bonds and 
negligible bond length fluctuations 

(3.14) 

~r -it\ f e'^e'^^ \ fa \ f e'^e'^" 
-it -r) \ e-'^e'^") \0 a* J \ e-'^e'^" 

V it\ fe'^ae'^ 
it r)\ e-'^a*e-'^ 



U{k) 



where 9i,...,9b are B independent random variables that are uni- 
formly distributed over the interval [0, 27r). The remaining parameters 
such as the central electron-electron scattering matrix a and the An- 
dreev scattering amplitudes tb are assumed to be flxed in one ensemble. 
We will call the ensemble homogeneous if all Andreev scattering am- 
plitudes are equal tb — t, otherwise we call it inhomogeneous. We 
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will now solve the homogeneous case and discuss some inhomogeneous 
ensembles numerically. 

3.2.1. The homogeneous ensemble. In this case the ensemble aver- 
aged spectral counting will turn out to be completely independent of 
the choice of the central electron-electron scattering matrix, and also 
of the number of bonds. In fact the case can be reduced to the case 
B = 1 that we have already solved. To see this we write 6b = d + S6b 
such that "^SOb = 0. We will first do the average over 9 and observe 
that this can be solved by diagonalisation. The straight forward result 
of the average over 6 will not depend on 66b. 
To do the average over 6 let us fix 50 and diagonalise 

(3.15) e^^Ve*''^ = Oe'"0^ 

where O is a real orthogonal matrix (remember that a is symmetric) 
and e**^ = diag(e*'^^, . . . ,6*^-^) contains the unimodular eigenvalues on 
the diagonal (we will omit the explicit dependence of the orthogonal 
matrix and the eigenvalues on 66b). We may now write 
(3.16) 

O [tt r [ g-i(2e-2«L+^+r.) I I 



U{k) 



and observe that the spectral counting function can be written as a sum 
over B modes each of which is equivalent to the single lead (the case 
B = 1). As a consequence the position and size of the spectral gaps 
are the same as for the single lead and the averaged density of states 
and the averaged spectral counting function are described by (12.121) 
and fl2.13p respectively (any dependence on the 66b has dropped). 
It is remarkable that this result is universal - it only depends on the 
Andreev scattering amplitude and on no other system parameter (the 
dependence on the length L is trivial). 

3.2.2. Some inhomogeneous ensembles. Different Andreev scatter- 
ing amplitudes at the peripheral vertices introduce an additional amount 
of disorder in the inhomogeneous case. One may expect that the form 
of the averaged density of states now depends on the actual choice of 
amplitudes tb and on the central electron-electron scattering matrix a 
- and the same applies to the existence and size of spectral gaps. In 
this section we will only discuss the extreme case tb G {0, 1} - that 
is, at some peripheral vertex one has either complete Andreev reflec- 
tions tfc = 1 or complete normal reflections t;, = 0. Let M < 5 be 
the number of bonds with total Andreev reflection such that t;, = for 
h = 1, . . . S - M and = I ioi h = B - M + I, . . . ,B. The mean 
probability to be Andreev scattered is = 
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One consequence of our choice of Andreev scattering amplitudes is that 
averaging over the phases 9b~m+i, ■ ■ ■ ,0b is completely irrelevant. To 
see this let us assume an electron is scattered into one of the bonds 
b = B — M +1, . . . , B. It accumulates a phase e^i^b+i^L) propagates 
from the centre to the peripheral vertex where it is scattered back as a 
hole and acquires a further phase factor —i. When the hole has propa- 
gated back to the centre the overall phase which has been accumulated 
is —ie^'^^ where the contribution from the random phase Oi, has been 
cancelled exactly. 

One further implication can be drawn from looking at the eigenvectors 
of the matrix a{0) — e'^cre'^ where d — diag(^i, . . . , 9b) is fixed. In 
general it is not possible to construct eigenstates of the Andreev star 
graph from the eigenvectors a because these get coupled in a differ- 
ent way at normal reflecting and at Andreev reflecting peripheral ver- 
tices. However, if d has an eigenvector a = (0, . . . , 0, os-m+i, ■ ■ ■ , os)'^ 
with eigenvalue e*" which is supported completely on the bonds h ~ 
B — M + 1, . . . , B then it is straight forward to show that U{0) has 

the two eigenstates ( _|_gia/2^ 1 which are also supported on the bonds 

b — B — M + 1, . . . , B. The corresponding eigenvalue is which 
contributes to the spectrum of the Andreev graph at = n + ^ (n 
integer). Furthermore, if there is such a state for some fixed choice of 
then there will be such a state for any element in the ensemble and one 
will observe a corresponding step in the averaged counting function. A 
similar statement docs not hold for states that are supported on the 
normal conducting bonds b = 1, . . . ,B — M. 

In the following we will restrict our discussion to two choices for the 
central electron-electron scattering matrix that are very well studied in 
the context of normal quantum graphs: 
i. the Neumann (or Kirchhoff) scattering matrix 

(3.17) CNeumannbb' — — ^bb' 

Jd 

and ii. the discrete Fourier transform (DFT) matrix 

(3.18) aj^FTbb' = ^e*'""^ ■ 

For quantum star graphs (with a scalar wave function and with incom- 
mensurable bond lengths) it is known that the DFT matrix generates 
universal spectral statistics as known from quantum chaotic systems 
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Figure 2. Dotted and full lines: averaged spectral 
counting function for Andreev graphs with negligible 
bond length fluctuations, inhomogeneous Andreev scat- 
tering, and a Neumann e-e scattering matrix. The two 
dotted lines are for B = 20, M = 6 (blue) and for 
5 = 80, M = 24 (green) such that pa = 0.3 in both 
cases. The full (black) line is the expected limiting func- 
tion and hardly differs from the line for B = 80. 
Dashed (red) line: averaged spectral counting function 
for the same inhomogeneous ensemble but with one set 
of randomly chosen bond lengths. 



and the Gaussian orthogonal ensemble (GOE) of random matrix the- 
ory. In contrast, the Neumann star graph belongs to a different uni- 
versality class with intermediate spectral statistics. 
We first discuss the inhomogeneous ensemble with a Neumann e-e 
scattering matrix crNeumann- The spectral decomposition of cxNeumann 
is straight forward. It has one eigenvalue 1 and B — 1 eigenvalues 
— 1. For our discussion it is important that one can construct M — 1 
orthogonal eigenvectors of ciNeumann that are supported on the bonds 
b = B — M + 1, . . . , B. According to our discussion above this leads 
to 2(M — 1) eigenvectors of the quantum map U{0) with eigenvalues 
±2 and the averaged spectral counting function will have steps of size 
= Pa - at = ±i -t- n (n = 0, ±1, ±2, . . . ). On each of these 
modes the probability to be Andreev scattered is unity. The remaining 
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2{B — M + 1) modes are hardly coupled (at least for large B) by An- 
dreev scattering - the mean probability to be Andreev scattered among 
these modes is of the order ^ and we may expect that the effect on 
these modes is negligible in the limit B ^ oo (where = ^ is fixed) 
such that the contribution to the averaged spectral counting function is 
linear. Figure [2] compares the numerically obtained averaged counting 
function for an Andreev graph with a Neumann e-e scattering matrix 
with the expected limiting function 

(3.19) A^(K)Neumann = (1 - Pa)k. h ^ Pa© ( " ^ 

n=— oo ^ 

Our expectation seems to be confirmed by the numerics and suggests 
that while there remains no gap in the averaged density of states the 
latter is reduced compared to Weyl's law near k = (modulo period- 
icity). 

Let us now look more closely at the spectrum near k = for a finite 
size of the graph. Figure [3] shows the integrated distribution /i(k) of 
the first positive eigenvalue of the Andreev graph as a function of ^^/t 
(that is we measure the spectrum in units of the mean level spacing 
according to Weyl's law). The numerics suggests that even on this 
blown-up scale there is no gap. /i(fi:) may either start algebraically 
Ji ~ hC or as a power series. In both cases it implies that there is a 
further reduction (in addition to the one discussed above) of the den- 
sity of states on this blown-up scale. We used 1000 realisations in our 
numerics which is not sufficient to resolve the region < 0.01 so 

our results are not fully conclusive. 

The behaviour of the spectral counting function and is very 

different if we choose the DFT e-e scattering matrix. Figure H] shows the 
numerically averaged spectral counting function for this case. There are 
several interesting features in this spectrum. Again, we find a jump at 
= I which is connected to eigenvectors of (Tdft which are supported 
on the bonds h = B~M+1, . . . ,B where the Andreev scattering occurs. 
However, the step is much smaller compared to the Neumann case and 
the probability for Andreev scattering on the remaining modes is large 
enough to have an effect. For large B the spectral counting function 
converges (numerically) to a limiting function - the plots in Figure H] 
are for B = 80 and hardly differ from the case B = 40 (not shown). 
Most importantly the counting function seems to drop to zero at a 
finite positive value of k. 

The numerically obtained distribution of the first eigenvalue /i(k) in 
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Figure 3. Integrated Distribution of the first positive 
eigenvalue of three Andreev star graphs with neghgible 
bond length fluctuations, inhomogeneous Andreev scat- 
tering {pa = 0.3 for all curves) and Neumann e-e scatter- 
ing matrix. The three plots correspond to an increasing 
number of bonds (B = 20, 40, 80). 

Figure [5] reveals that a spectral gap appears at least in the limit B oo 
(at constant pa)- 

3.3. Andreev star graphs with random bond lengths. Let 

us now discuss the influence of bond lengths and discuss ensembles 
of Andreev graphs for which the bond length fluctuations cannot be 
ignored. For a graph with B bonds let us choose a fixed set of bond 
lengths in the interval < L;, < 1 and now average over the Fermi 
wave number kp- The main effect of the different bond lengths is 
that one cannot get the spectrum of the graph from just diagonalising 
the quantum map at k = and the spectral properties are no longer 
periodic. In fact one should expect that any features in the counting 
function which differ from Weyl's law can only survive near k = 0. 
We have calculated the averaged counting function numerically using 
the trace formula fl3.1ip . We chose the bond lengths independently in 
the interval < Lf, < 1 with a random number generator. All plots we 
show are for a single realisation of the bond lengths (we have assured 
that the figures are generic by comparing to other choices). Indeed 
Figures [2] and H] show that the averaged counting function for both 



16 HOLGER FLECHSIG AND SVEN GNUTZMANN 



1.0 

{N{k)) 

0.75 



0.5 



0.25 








0.2 0.4 0.6 0.8 1.0 



Figure 4. Full (blue) lines: averaged spectral counting 
function for Andreev graphs with negligible bond length 
fluctuations, inhomogeneous Andreev scattering, and a 
DFT e-e scattering matrix. The fat line corresponds to 
5 = 80, M = 24 {pA = 0.3) and the thin line to B = 80, 
M = 40 {pa = 0.5). Plots for i? = 40 and the same 
values for pa (not shown) are almost identical. 
Dashed (red) lines: averaged spectral counting function 
for the same inhomogeneous ensemble but with one set 
of randomly chosen bond lengths. 

the Neumann and the DFT e-e scattering matrix has lost all special 
features when — k is on the order of one. For smaller values — k <C 1 
the fluctuations in the bond lengths are not effective and one gets back 
to the case discussed before. In short, the bond length fluctuations 
destroy some special features of the previous ensembles with negligible 
bond length fluctuations. The bond length fluctuations do not have a 
strong influence on the reduction of the density of states near k = or 
the existence and size of a gap, at least when ^ 1 @- 

4. Conclusions 

We have shown that the spectra of Andreev star graphs have a 
variety of special spectral properties that are related to the introduction 
of Andreev scattering into a quantum graph. There is a spectral gap in 



This should be taken with care since in our model the bond length fluctuation 
are still moderate. 
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Figure 5. Integrated Distribution of the first positive 
eigenvalue of Andreev star graphs with neghgible bond 
length fluctuations, inhomogeneous Andreev scattering 
and DFT e-e scattering matrix. The three plots on the 
left correspond to pa = 0.3, the three on the right to 
Pa = 0.5. In both cases we have plotted the distribution 
for B = 40 (dotted red lines), B = SO (dashed green 
lines), and B = 160 (full blue lines). 

homogeneous ensembles of Andreev graphs. Both the spectral counting 
function and the spectral gap are universal and only depend on the 
probability of Andreev scattering pA- 

Inhomogeneous ensembles contain additional disorder and the universal 
results from the homogeneous case are not applicable. Instead there is a 
rich variety of phenomena that depend on the properties of the central 
scattering matrix. For some special choices we could show numerically 
that the spectral counting functions converge to a limiting function as 



This first numerical investigation indicates that the appearance and 
size of a spectral gap is strongly related to the localisation properties 
of the eigenstates of the central e-e scattering matrix. 
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